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ABSTRACT 

We present a detailed Fermi/LAT data analysis for the broad-line radio galaxy 3C 120. 
This source has recently entered into a state of increased 7-ray activity which mani¬ 
fested itself in two major flares detected by Fermi/LAT in September 2014 and April 
2015 with no significant flux changes reported in other wavelengths. We analyse avail¬ 
able data focusing our attention on aforementioned outbursts. We find very fast vari¬ 
ability timescale during flares (of the order of hours) together with a significant 7-ray 
flux increase. We show that the ^ 6.8 years averaged 7-ray emission of 3C 120 is likely 
a sum of the external radiation Compton and the synchrotron self-Compton radiative 
components. To address the problem of “orphan” 7-ray flares and fast variability we 
model the jet radiation dividing the jet structure into two components: the wide and 
relatively slow outer layer and the fast, narrow spine. We show that with the addition 
of the fast spine occasionally bent towards the observer we are able to explain observed 
spectral energy distribution of 3C 120 during flares with the Compton up-scattered 
broad-line region and dusty torus photons as main 7-rays emission mechanism. 

Key words: galaxies: individual (3C 120) - galaxies: jets, general - radiation mech¬ 
anisms: non-thermal - gamma-rays: galaxies 


1 INTRODUCTION 

The recent detection of 7-rays from several broad-line ra¬ 
dio galaxies (BLRGs: 3 C 120, 3C 111, s-iid Pictor A) 
by F e rmi/LAT telescope (lAbdo et al.l l20ld : iKataoka et al.l 
I2OIII : iBrown fc Adama I2OI2I I constituted this sub-class of 
active galactic nuclei (AGN) as high energy emitters. To- 
gether with radio galaxies such as M87 (lAbdo et al.ll200^ 1 or 
Gen A llAbdo et al.ll201ol l seen by Fermi/LAT earlier they 
form an interesting group of AGN where high energy ra¬ 
diation is observed despite the absence of strong Doppler 
boosting as we observe in typical blazars. These non-blazar 
AGN (dubbed “misaligned blazars”), although much fainter 
due to geometrical effects, are very interesting objects as 
they reveal both accretion disk radiation (at least in case of 
BLRGs) and jet emission in their broadband spectra. 

From infrared to hard X-rays BLRGs show typical 
thermal emission related to accretion disk around super- 
massive black hole, which can be further divided into di¬ 
rect accretion disk radiation peaking in optical wavelengths, 
dusty torus infrared emission and power law non-thermal X- 
ray component most probably originating from disk corona 
dZdziarski fc Grandill200ll . iGrandi fc Palumbo! lioOTl f. Non- 
thermal radio and 7-ray emission is thought to originate 
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from relativistic jet and is being only weakly Doppler 
boosted due to large viewing angles Sobs 10°. 

3C 120 is a nearby {z = 0.0 33) BLRG wi th Fanaroff- 
Riley type I radio morphology dWalker et al.l I2OI5I I . The 
source has been actively monitored in all wavebands from 
radio up to X-ra ys and recent detection by Fermi/LAT 
llAbdo et al.ll201(J 'l made it possible to study accretion disk 
and jet interaction by modelling its broadband spectra. 

In this article we present detailed data analysis of high 
energy data for the whole Fermi/LAT dataset (~ 6.8 years) 
and we especially focus our attention on two major outbursts 
that happened in September 2014 and April 2015. We study 
both temporal and spectral properties of these flares. 

This paper is organised as follows. In section [2] we 
present the details and results of Fermi/LAT data analysis. 
In section[3]we present discussion concerning the location of 
the 7-ray emission region (the “blazar zone”), the adopted 
“spine-layer” jet model as well as details on spectral mod¬ 
elling. Summary and conclusions are presented in section [I] 


2 FERMI/LAT DATA ANALYSIS 

The 3G 120 7-ray spectra and light curves within 100 MeV 
to 100 GeV energy range were obtained by analysing 
about 6.8 years of Fermi/LAT data from August 4th, 2008 
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to May 26th, 2015 (hereafter; all Fermi/LAT data). We 
reduced the data using Fermi Science Tools package ver¬ 
sion v9r33p0 selecting only source class 2 events and us¬ 
ing P7REP_S0URCE_V15 instrument response function. Given 
a rather soft spectrum of 3C 120 and a large instru¬ 
ment point spread function (3.5° at 100 MeV) we selected 
events within 20° (ROI) from nominal source radio posi¬ 
tion (R.A.=68.296°, Dec.=5.354°) to properly model the 
contribution of other sources. To avoid contamination form 
Earth’s limb photons we also applied a zenith angle > 100° 
cut. 


We used the standard gtlike tool to mod el the ROI 
using maximum likelihood method dCas ailizi). Model in¬ 
cludes 69 po int sources from LAT 4 -year Point Source Cat¬ 
alog l3FGL. rAckermann et al1l2015r) as well as Galactic and 
isotropic diffusion emission templates Q. We model both 
flux normalisation and power law spectral index for sources 
within a radius of 10° from the ROI centre. Other sources 
within the ROI are fixed to their 3FGL values. After initial 
fit, all sources with test statistics (TS) lower than 1.0 were 
excluded from the model and the procedure was repeated 
until convergence. 

Two sources in the 3FGL are located fairly close to 
3G 120 i.e. 3FGL J0432.5-b0539 and 3FGL J0426.6-b0459. 


The latter is 1.66° away yet the former is only 0.35° away 
from the source of interest having photon index of a = 2.7 
which makes distinction between the two so urces unclear. 
Howev er, by fitting 7 -ray 3G 120 position iTanaka et al.l 
(l2015ll found that it lies 0.028° ± 0.088° from its radio po¬ 
sition. Also, 3FGL J0432.5-1-0539 LAT position was deter¬ 
mined with a 95% error of 0.15°. Therefore, we include all 
three sources in the model treating them as separate. 

For the whole dataset (all data) we used binned 
likelihood method. 3G 120 was detected with TS = 
156, a 100 MeV - 100 GeV flux of E = (3.0 ± 0.6) x 
10 “® ph cm“^ s~^ and a power law photon index a — 
2.71 ± 0.09. We calculated light curve and spectra by di¬ 
viding the data into time bins (30-days, 83 bins) or energy 
bins (6 equal bins in logarithmic scale covering the chosen 
energy range), and by applying previous procedure to model 
the ROI with both flux normalisation and power law index 
as free parameters. In case the TS value was lower than 9 
(corresponding to significance cr ~ 3) 95% upper limits were 
calculated assuming power law index of a = 2.7. Figure [T] 
presents source 7 -ray spectra. 

Figure [2] presents 3G 120 light curve. Source was ob¬ 
served with significance larger that 3cr several times in 
recent years. We focus our attention on two recent out- 
bursts: in Septemb er 2014 (hereafter: flare 1; reported by 
iTanaka et al.f 201^) and in April 2015 (here after: flare 2; 
reported by iFermi/LAT CollaborationI l2015l L To analyse 
Fermi/hAF data during flare 1 and 2 periods we used un¬ 
binned likelihood analysis method to properly account for 
low photon count numbers. During flare 1 in September 2014 
(30-days of data) 3G 120 was observed with TS = 11, power 
law photon index a = 2.50 ± 0.32 and its 95% flux upper 
limit was estimated as Ful = 1.2 x 10“^ ph cm“^ s“^. Fig¬ 
ure [ 3 ] (left panel) presents light curve for September 2014 
with clear flux increase on 24th September. Right panel of 


gll_ieni_v05_revl. f it and iso_source_v05.txt, respectively. 


the same figure presents light curve for that day with 2-hour 
time bins. On 24th September 3C 120 was detected with 
TS = 54 with photon flux F = (9.2T2.7) x 10“^ ph cm“^ s~^ 
and power law index of tr = 2.34 ±0.25. Even for the 2-hour 
highest flux time bin TS = 28 indicating detection above 5a 
significance level. 

We applied the same procedure to investigate flare 2. 
Figure [ 4 ] presents monthly (1-day time bin) and daily (30- 
minutes time bin) light curves for April 2015 and 24th April 
2015, respectively. Monthly averaged flux is F = (1.45 ± 
0.09) X 10~^ ph cm“^ s“^ with power law index a = 2.47 ± 
0.03 (TS = 57) and daily averaged flux for 24th April is 
F = (2.6 ± 0.35) X 10“® ph cm“^ s“^ with power law index 
a = 2.22T0.11 (TS = 281). Source was also detected within 
30-minutes time bin with large TS = 89. 

Figure [ 5 ] presents flares 1 and 2 7 -ray spectra averaged 
over outburst day and shortest time bin with detection above 
3cr i.e. 2 hours for flare 1 and 30 minutes for flare 2. 

For flaring periods we additionally checked whether 
they could be associated with nearby sources 3FGL 
J0432.5±0539 and 3FGL J0426.6±0459. Light curve anal¬ 
ysis for these sources indicated no correlation between their 
flux increase and presented 3C 120 outbursts. 

Table [T] summarises Fermi/LAT data analysis results. 


3 DISCUSSION 


3.1 Location of the blazar zone 

FermilhAJ: data suggest that 3C 120 underwent very rapid 
changes in 7 -ray flux. Both aforementioned flares happened 
on very short timescales. Assuming that characteristic vari¬ 
ability timescale Uar equals the shortest time bin for which 
the source has been detected with significance larger than 
3a we estimate fvar = 2h for flare 1 and tvar = 30min for 
flare 2. In a similar way 3C 120 7 -ray variability timescale for 
six years of Ferm i /hAT data was estimated to be 5 —10 days 
by ITanaka et al.l (120151 1. 

Variability timescale and causality arguments allow us 
to estimate radius of the emission region R ^ cFUar where 
V = [r(l — /3cos 0obs)l~^ is the Doppler factor, F is the jet 
Lorentz factor, /3 is the jet velocity in the speed of light units 
and 9obs is the angle towards th e observer. For 3G 120 the 
average value of F = 5.3 ± 1.2 (Ijorstad et al.l 1200^ 1 which 
is in agreement with further work by Casadio et al.l ll2015h . 
The angle towards the observer Oobs was estimated in several 
papers using V LBI observations to have differ ent mean val- 
ues: from 9.7° dHovatta et al.|[2009l l to 20.5° (Ijorstad et al.l 
I 2 OO 5 II indicating that the jet might be changing its direc¬ 
tions or that the radio knots might be moving not along the 
jet axis. For further considerations we adopt a mean value of 
^obs = 15°. Assuming aforementioned variability timescales 
and the conical jet opening angle djet = 1/r we estimate the 
radius of the emitting region to be F < 1.0 x 10^^ (11/5.3) cm 
for all data, R < 1.1 x 10^®(Il/5.3) cm for flare 1 and 
F < 2.9 X 10^'‘(F>/5.3) cm for flare 2. 

Using F = 5.3 and Oobs = 15° which results in 
F = 3.7, in agreement with values reported in the litera- 
ture (iHovatta et al.l[2009l: Ijorstad et ^l2005l: ICasadio et al.l 

l2015ll . we estimate the location of 7 -ray emission region 
(the blazar zone) from the central black hole (BH) at 
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Figure 1. 3C 120 spectra in 100 MeV — 100 GeV energy range for the whole dataset from 4th August 2008 to 26th May 2015. Arrows 
indicate 95% upper limits for detection with significance lower than 3 (t. 




Figure 2. 3C 120 light curve (100 MeV—100 GeV) from 4th August 2008 to 26th May 2015 with 30-days time bins. Red points correspond 
to detection above 3 (t significance while grey points indicate 95% upper limits for detection with significance lower than 3cr. Blue, dotted 
line and light-blue shaded area correspond to 6.8 years average FermifLAT flux and flux error, respectively. 


r — FR < 5.3 X 10^^ cm « 0.17 pc for all Fermi/LAT 
data. Similar es t imate s were calculated by other works. 
ICliatteriee et al] ll2009l) calculated a de-projected distance 
of ~ 0.5 pc from the central BH to the 43 GHz VLBA 
radio core by using the time lag between the dip in the 
X-ray flux (assumed to originate in accretion disc corona) 
and the ejection of superluminal component from the radio 
core. By analysing Fermi/LAT lighcurves and radio data 


ICasadio et al.l ll2015l) found that 7 -ray emission region lies 
about ~ 0.13 pc upstrea m from the rad i o core . Similar con¬ 
clusions were drawn by I Tanaka et al.l ll2015l) locating the 
emission region at 0.1 — 0.3 pc from the central BH. 

By assuming that the 7 -ray emission region is transver- 
sly uniform relation r = FR gives the location of the blazar 
zone at r < 5.8 x 10^® (11/5.3) (r/5. 3) cm for flare 1 and 
r < 1.5 X 10^® (11/5.3) (r/5. 3) cm for flare 2. 
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Figure 3. 3C 120 light curve (100MeV — 100GeV). Left panel: data for September 2014 with 1-day time bins. Right panel: data for 
24th September 2014 with 2-hour time bins. Red points correspond to detection above 3 (t significance while grey points indicate 95% 
upper limits for detection with significance lower than 3(7. Blue, dotted line correspond to monthly (left panel) and daily (right panel) 
average FermilhAT flux. 
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Figure 4. 3C 120 light curve (100 MeV — 100 GeV). Left panel: data for April, 2015 with 1-day time bins. Right panel: data for 24th 
April 2015 with 30-minutes time bins. Red points correspond to detection above 3(7 significance while grey points indicate 95% upper 
limits for detection with significance lower than 3(7. Blue, dotted line correspond to monthly (left panel) and daily (right panel) average 
FermilLAT flux. 


3.2 Spine-sheath geometry 

Outbursts observed in 3C 120 by Fermi/LAT are extreme 
both in terms of temporal and spectral characteristics. Both 
flares happened at very short timescales, of the order of 
hours, which is unusual as no such rapid outbursts have been 
observed before in 7 -rays for any source of this kind. Intrest- 
ingly, there are no reports on increased activity for other 
wavelengths. For flar e 1 several f ollow-up obs ervations were 
reported: in optical llNescil[2014h . in X-rays (iLohfink et al.l 


l2014h and in near infrared (iHasan et al. IliH however none 
of them reported any unusual source behaviour. All of them 
were carried out days after flaring activity hence it is difficult 
to directly connect those observations to 7 -ray brightening. 
No significant change in flux ha s been observed in so ft X- 
rays with the MAXI telescop(j3 (iMatsuoka et al.|l2009l) and 
in hard X-rays with the BAT instrument on board the Swift 

^ http://maxi.riken.jp/top/index.php?cid=l&jname=J0433+053 
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Figure 5. 3C 120 spectra in 100 MeV — lOGeV energy range for flare 1 (left panel) and flare 2 (right panel). Arrows indicate 95% upper 
limits for detection with significance lower than 3 (t. Red points present daily averaged spectra while grey points show spectra for highest 
flux 2-hours time bin (flare 1) and 30-minutes time bin (flare 2). 


Table 1. Summary of Fermi/LAT data analysis results for 3C 120. a and TS denote power law photon index and test significance, 
respectively. 


Time period 

Flux [10 ^ph cm ^ s 

a 

TS 

6.8 years 

0.3 ±0.06 

2.71 ± 0.09 

156 

Flare 1 

September 2014 

1.2 (95% upper limit) 

2.50 ±0.32 

11 

24th September 2014 

9.2 ±2.7 

2.34 ±0.25 

54 

2-hour highest flux time bin 

21 ±7.5 

2.40 ±0.35 

28 

Flare 2 

April 2015 

1.45 ±0.09 

2.47 ±0.03 

57 

24th April 2015 

26 ±3.5 

2.22 ± 0.11 

281 

30-minutes highest flux time bin 

58± 14 

2.07 ±0.21 

89 


satellitflEI l|Barthelmv et al]l2005h . for both flaring periods. 
We have no data on flaring activity in other wavelengths 
coinciding with flares 1 and 2 periods. 

3C 120 7 -ray flux increase on very short timescale, ob¬ 
served spectral hardening, no flaring behaviour in X-rays 
and no indication of rapid changes in other wavelengths 
constitute a challenge for inter pretation of those facts in 
terms of usual disk-|-jet models (iTanaka et al.ll^l5^ . No in¬ 
dication of increased activity except for y-rays suggest that 
high energy component and radiation at other wavelengths 
may be produced at different locations. This observation 
also supports the idea that BLRGs broadband spectra con¬ 
sist of disk-related and jet-related components. However, no 
change in hard X-rays during 7 -ray flaring sets certain limi¬ 
tations on radiation mechanism responsible for high energy 
component in relativistic jet i.e. assuming one-zone model 
for production of radiation in the jet, 7 -ray emission is con¬ 
strained by the constant X-ray flux. 

® http://swift.gsfc.nasa.gov/results/transients/weak/ 
3C120.1c.txt 


7 -ray flux during flare 1 was ~ 50 times and during 
flare 2 ~ 150 about times larger than the average flux for 
3C 120 for all Fermi/LAT data. Total apparent luminosity 
from the jet is L = (T>^/r)77rad7e»7dissFjet where ? 7 diss is the 
total jet power dissipation efficiency, 7 e describes what part 
of dissipated energy is transferred to electrons, ? 7 rad is the ra¬ 
diative efficiency and Ljet is a total jet power before dissipa¬ 
tion. Increase of total apparent luminosity can be achieved 
by either increase in jet bulk Lorentz factor T, more effi¬ 
cient jet power dissipation/radiative energy dissipation or 
increase in the total jet power. Let us focus on possible im¬ 
plications of either of these possibilities by assuming that 
the low-energy jet spectral component is synchrotron radia¬ 
tion whereas high-energy jet spectral component is external 
radiation Compton emission (ERC) (see section [3.311 . The 
ERC luminosity can be estimated as Lerc ryLjetT^ by its 
dependence on the jet power (where rj = 77rad?7e??diss) but 
also as Lerc oc oc Ldisk where is the external radi¬ 
ation energy density in the jet co moving frame which is a 
function of the accretion disk luminosity Ldisk. Synchrotron 
radiation luminosity can be estimated in a similar manner as 
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Lsyn oc Mb oc Lb oc crLjet/(l + O') where Mb is the magnetic 
field energy density in the jet frame, Lb is the magnetic 
field energy flux and a is the jet magnetisation defined as 
a — Lb /Lkin where Lkin is the jet kinetic energy assumed to 
be dominated by cold protons. In a one-zone leptonic model 
LERc/Layn OC Ue^tT^/u'g thus we have 


1 -I- a 



V- 


If the increase in 7 -ray band ALbrc was caused by 
the increase in the jet power ALjet this would lead to 
even more significant increase in synchrotron luminos¬ 
ity i.e. ALsyn oc (ALjet)^. Increasing overall energy dissi¬ 
pation efficiency i.e. ALbrc oc Arj results in ALgyn oc 77. 
Finally, as the expression above is not a function of F, an 
increase of the F factor causes significantly larger 7 -ray flux 
[ALbrc oc (AF)^] while it does not affect synchrotron lumi¬ 
nosity. Such simple analysis proves that first two possibilities 
have to be excluded as they would lead to increase in not 
only 7 -rays but also other wavelengths which is not observed 
in case of 3C 120 flaring states. Third possibility can also 
be discarded because increase in F factor would lead to nar¬ 
rower Doppler cone lowering total apparent luminosity at all 
wavelengths. A straightforward conclusion is that to explain 
such “orphan” flares in 3C 120 there is a need for a sepa¬ 
rate emission zone. Such a zone should satisfy the following 
conditions: 


• to have the ability to occasionally produce large 7 -ray 
fluxes, strongly beamed towards the observer, 

• to have very large F Lorentz factor to boost ERC radi¬ 
ation avoiding too large synchrotron luminosity, 

• be compact enough to address the very fast variability 
during flares. 


We propose that such a component might be represented 
by the fast moving and wiggling spine. Spine-sheath/layer 
jet structures have been considered in the literature in va¬ 
riety of models propos ed to explain specific properties of 


jetted objects (see e.g. Celotti et al.l 200^ ^fiiseUiniet^^y 


2 OO 5 I: iTavecchio fc Ghisellinil 20081 : D’Arcangelo et ^ 20091 : 


Mimica et al. IM) . In the case of a magnetically-arrested 


disk (MAD) scenario of AGN jets launch ing, which is likely 
to apply to 3C 120 (iLohfink et al.l[2013lf . faster spine and 
slower layer can naturally result from non-unifo rm mass 
loading of a jet at its base (iMcKinnev et al1l2012| j. In such 
a case mass loading is driven by interchange instabilities. 
Because these instabilities are non-axisymmetric, they may 
also be responsible for wiggling of a jet and its velocity varia¬ 
tions. Jet dir ection can be also altered by t he current-driven 
instabilities (iNalewaiko fc Begelmanl|2012l and refs, therein). 
In this work we follow the spine-layer idea to model the spec¬ 
tra of flares 1 and 2 in 3C 120. 

We assume that the conical jet with opening angle 6jet is 
stratified in such a way that it is composed of two elements: a 
fast conical spine with bulk Lorentz factor F’®, opening angle 
^iet and a conical layer (“sheath”) with bulk Lorentz factor 
F* = F and opening angle ^jet — ^jet- The layer forms the 
main jet body being aligned with its axis and the spine forms 
an addition inside the jet. We assume that F® 3> F* and 
^let < ^]ef We also assume that, while layer stays straight, 
the spine is subjected to change its direction with respect 
to jet axis and, therefore, with respect to the observer. In 


other words we assume layer’s Doppler factor is constant 
while spine’s Doppler factor is variable due to variable 
angle towards the observer ^obs- We associate this variability 
timescale with observer 7 -ray variability. 

Note that we do not discuss the origin of spine-layer 
configuration nor do we calculate spine-layer geometry from 
basic principles e.g. using jet formation theories. However, 
since we assume that the layer actually forms the underlying 
and stable base of the jet its parameters, 0]^^ and F* can 
be inferred directly from observations. Spine’s parameters 
on the other hand are model free parameters i.e. they are 
adjusted to model source spectra. 

Figure [ 6 ] presents a schematic view of the geometry of 
the presented model. 


3.3 7 -ray emission mechanisms 


While the low energy jet radiation component is likely syn¬ 
chrotron emission the radiation mechanism for 7 -rays re¬ 
mains unclear. P ossible mechanisms inc l ude self-synchrotron 
Compton (SSC; iMaraschi et al.l 1 19921 : iBloom fc Marscto 

S or external radiat ion Compton lERC: IPermer et al. 

ISikora et al.l Il994l l on broad-line region (BLR) and 
hot dust region (HDR) seed photons. 

IPozo Nunez et IdJ (l2014l'l calculated the location of the 
thin-disk BLR in 3C 120 between 22 and 28 light days (l.d.) 
from the central BH with reverbe ration mapping tec h nique . 
Similar results very obtained by iKollatschnv et aP (l2014l l 
who found BLR to be stratified with He H emission line lo¬ 
cated at 12±7 l.d. and Ha line much further at 28.5±8.5 l.d. 
away from the BH. For the purpose of this work we will 
adopt a mean BLR radius value of 25 l.d. which corresponds 
to rsLR = 6.5 X 10^® cm ~ 0.02 pc. Recent models of duty 
torus and BLR in AG N (IC zernv fc Hrynie wicd I20l11j and 
IR torus observations (|Kishimoto et al.ll201l| j suggest that 
the HDR characteristic radius thdr ~ lOrBLR thus we set 
rHDR = 6.5 X 10^^ cm « 0.2 pc. External photon fields en¬ 
ergy density in the jet co-moving frame for r larger than 
the characteristic radius is Mext = ^CcFF^Ldisk/ 47 rr^c where 
fcF is BLR (HDR ) covering factor (we assume ^cf^ = 0-1 
I Sikora et al.il200^ and ^cf^ = 0.3 due to larger torus size) 
and ^ « 0.1 is a facto r accounting for fla t geometry of pho¬ 
ton emission regions ll Janiak et al .11201^ . 3C 120 accretion 
disk flux at ~ 10 eV is ~ 1.5 x 10“^® ergcm^ s“^ result¬ 
ing in disk luminosity Laisk ~ 1.7 x 10^® e rgs~^ after ap¬ 
plying bolometric correction of ~ 4.5 from [Richards et al.l 
( 20061). Similar v alue of disk luminosity was obtained by 
lOgle et 'aD (l2005l L We note that this value and our assump¬ 
tion concerning dusty torus geometry is consistent with ob¬ 
servation that the inner edge of dusty torus in AGN is lo¬ 
cated approximately at graphite sublimation radius r sub = 
1.6 X ~ 6.6 X 10^^ cm (iMor fc Netzer|[201^ . 

To find an efficient and dominant mechanism for 7 -ray 
production in case of all data averaged spectra let us esti¬ 
mate the ratio between ERC and SSC luminosities produced 
in a uniform jet. For aforementioned values we obtain 


^ext 


3 X 10“® 



-) 7 — 

5.3 y V’-hdr/ 


erg cm 


Synchrotron radiation ene rgy density in the jet frame is 
Wsyn = Lsyn/'^'!^'D'^^fetr^c (Ijaujak et al.l l20f5l ’l - with syn¬ 
chrotron bolometric luminosity Lsyn ~ 1.6 x 10 ^^ erg s“^ 
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(see section EH and Fig. [2l) we get 


u 


The ratio between 
mately 


1.5 X 



-4 



erg cm 

.rnoR/ 

ERC and SSC luminosity is approxi- 


Tbrc 

Lssc 



(0jet/O.19)' , 


therefore ERC mechanism is dominating high energy radia¬ 
tive output. However, we note that this result strongly de¬ 
pends on Doppler factor D and jet opening angle 0jet - both 
these factor are poorly determined. Lowering the Doppler 
factor to a value of 2.4 (see section EU or lowering 9jet to 
a value of ~ 0.02 may result in indi c ation that SSC mech¬ 
anism is favoured [see I Tanaka et aP (l2015l ')1. Thus, despite 
our estimation we conclude that both ERC and SSC mech¬ 
anisms can be equally viable to explain the averaged data 
for 3C 120. 

Our conclusion has to be altered in the case of observed 
7 -ray flares where we use the sum of the spine and layer ra¬ 
diation to model source spectra. In such case production of 
7 -ray emission in spine via SSC mechanism is very unlikely 
as it would impose significant flux increase in radio (due 
to synchrotron emission) and X-ray band (due to spectral 
broadness of SSC component). ERC mechanism on external 
photons from BLR(HDR) is therefore favoured as it can ac¬ 
count for increased 7 -ray emission without the increase of 


flux in X-rays (see 3.2) provided its production in separate 
component (spine). 

In the spine-layer geometry model, during the period 
of source quiescence the total radiative output is dominated 
by the layer emission. Indeed, relatively low T factor and 
large viewing angles cause the jet emission to be significantly 
lower than the accretion-related emission. At this period of 
time the inner, fast spine may or may not exist - this fact 
is not resembled in the source spectra due to even larger 
viewing angles and spine power being lower than power car¬ 
ried by the layer (see 3.2). However, whenever the spine 
bends towards the observer (see Fig. |6] and section 3.2) the 
spine emission starts to dominate the total radiative out¬ 
put. In terms of spectral characteristics the ratio between 
ERC (from spine) and synchrotron (from layer) emission is 
-t'ERc/^ayn « «?xt/“B « X Wext/uR i-S- T is largely dom¬ 
inated by spine ERC radiation due to its much larger F® 
factor. 

Our estimates presented in previous paragraph (assum¬ 
ing values of Doppler factor observed for 3C 120) locate the 
blazar zone a t sub-parsec s cales deep inside BLR. In fact, 
as shown by Ijaniak et al.l ll2015h at distances lower than 
5 X 10^® cm from the central BH jet environment is domi¬ 
nated by photons from the accretion disk. In such circum¬ 
stances 7 -ray absorptio n via pair production is in evitable al¬ 
ready at GeV energies (iPoutanen fc Sternll201(]tl . The bent 
spine scenario, however, may provide a much larger size of 
the estimated emission region because of larger values of 
both Doppler factor H and Lorentz factor F (see 5.1) i.e. the 
estimated size R can be increased maximally by a factor 
(r®/r*)^ (assuming 6obs = l/F®) with respect to the value 
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obtained with the assumption that the emission region is 
calculated with full conical jet scenario. 

The spine-layer geometry essentially introduces a sec¬ 
ond emission zone and, therefore, inevitable interaction be¬ 
tween radiation produced in both spine and layer regions. 
The most significant radiative component to take into ac¬ 
count is layer’s (spine’s) synchrotron radiation Compton 
up-scattered in the spine (layer) region. However, due to 
geometry synchrotron radiation produced in spine would be 
significantly diluted in the much larger layer emission region 
dChisellini et al.l2005h . Therefore, it is most important to es¬ 
timate the spine Compton radiation on layer’s synchrotron 
photons. The ratio between ERC emission in spine Terc 
and external synchrotron Compton emission in spine Tesc 
is 

L s ^ /s 7^82^ 

ERC ^^ext i Wext 

- CK —;- OC =-, 

-^ESC ^syn.l 

where Usyn,i denoted layer’s synchrotron radiation energy 
density in spine and T”* = r®r*(l — /3° /3*) is an effec¬ 
tive L orentz factor between spine and layer dChisellini et al.l 
[200i). It can be shown that P' « r72r' for T" > T*. Using 
values calculated before we get Terc/^esc ^ '^0 hence spine 
Compton up-scattered synchrotron radiation from layer is 
much less luminous as compared to spine’s ERC radiation 
and may be omitted in further considerations. 


3.4 Model parameters and spectral modelling 


For spect ral modelling we us e numerical code already de¬ 
scribed in I Janiak et al.l ll2015ll . We point out main assump¬ 
tions of the numerical model below. 

Jet energy dissipation takes place in an active region of 
conical jet, at distances between [r, 2r] from the central BH 
and production of radiation is followed further up to 20 r. 
Electron evolution and pr oduction of radiatio n is calculated 
in a steady-state manner (ISikora et al.ll20f^ . 

We follow the evolution of electron energy distribution 
in the region of intere st by solving the kine tic equation for 
relativistic electrons (iModerski et al.l l2003l l which can be 
presented in the form 


dN^jr) 

dr 


A 

d7 



-b 


Q-ijr) 

cPV 


where N~i is the number of electrons per energy bin, /3 = 
— 1/r, d 7 /dr = (d7/dt')/(/3cr), d 7 /dt' are the elec¬ 
tron energy loss rates as measured in the jet co-moving 
frame, and Q-f{r) is the broken power-law electron injection 
function defined within [7inin, 7max] energy range, with spec¬ 
tral indices pi and p 2 . Injection function normalisation and 
electron spectrum break energy 7 b are calculated based on 
several parameters i.e. jet magnetisation a, jet pair content 
riejnp, jet production efficiency rjj, fraction of energy trans¬ 
ferred to electrons rjs, jet magnetic field B and energy dis¬ 
sipation efficiency » 7 diss. These parameters also allow to cal¬ 
culate jet energetics i.e. total jet power Tjet, pr otons kinetic 
energ y Lp and magnetic energy flux Lb (see I Janiak et al.l 
120151 for details). 

Synchrotron, SSC and adiabatic electron energy loss 
rates are calcu l ated using the procedure presented by 
IModerski et al.l (l2003l 'l. ERC energy loss rates and lumi¬ 
nosities of all radiation mechanisms as well the descrip¬ 


tion of the adopted here planar model of external radiation 
sources BLR and HDR were described in our previous works 
(ISikora et al.ll 20 f^ . I Janiak et al.ll20i5h . 

For all FermifLAT data we set the location of the blazar 
zone at r = 2.6 x 10^^ cm (see 3.1) assuming T = 5.3 
and Sjet = 1/ r. We assume the centr al BH mass Mbh = 
4.6 X 10® Mq iIPozo Nunez et ^|2014^ which for accretion 


cretion disk radiative efficiency 77disk = 0.1 results in ac¬ 
cretion power Mc^ ~ 1.7 x 10"*® ergs“^ and Eddington ra¬ 
dio of ~ 0.03. We note that BH mass estimation is very 
uncertain in 3C 120 as many works repor ted much lower 
values: 5.5 x IOTUq JPeterson et al.l[2Q )4 ), 3.0 x lO’^ 
( [hfarshall et al.l 120091 1 or 5.7 x lO’^M q ( Pozo Nunez et S 
I 2 OI 2 II . To model accretion-related emissi on in 3C 120 we use 
a radio-loud quasar radiation template dShang et al.ll2ojTh 
with jet radio emission truncated at ~ 10^® Hz. We match 
the template with th e X-ray data in 1 — 6 keV energy band 
jLohfink et al.l^2014^ . 


To estimate other input parameters of our jet emission 
model we choose them so that the calculated spectra match 
the observations. Having set an average value of the angle 
towards the observer 0 obs = 0.26 ~ 15 ° normalisation of 
calculated emission depends on total jet radiative efficiency, 
fraction of jet power effectively channelled to electrons in 
dissipation region and initial power carried by jet itself. Ra¬ 
diative efficiency at sub-parsec scales i.e. inside dense pho- 
ton fields from BL R and HDR was estimated to be > 0.5 
ll Janiak et ahllioisl l and initial jet power cannot significantly 
exceed accretion power. As 3C 120 is an FR I type of ra¬ 
dio source energy dissipation rjdms cannot be too low. The 
upper limit on rjdiaa can be inferred from the fact that the 
spectral peak zzpeak in 7 -rays is located at energies lower than 
100 MeV (see Fig. [TJ. For fixed values of injected electrons 
spectral indices the break in broken power law spectrum 
7 break OC 7 diss- SiuCe /Zpeak ^ U 7 break^ext We find that the 

upper limit on 7 break is ~ 10® (assuming ERC process on 
BLR photons with i/ext = J^blr = 10 eV). We assume that 
energy dissipation efficiency is 7 diss = 0.3 (which results in 
7break ~ 6 X 10^) aud half of that energy is transferred to 
electrons. 

For injected electron energy spectrum we set spectral in¬ 
dices to pi = 0.5 for 7 ^ 7 break aud P2 = 2.4 for 7 > 7 break. 
We set the minimum electron energy to 7 min = 1.0 and the 
maximum electron energy to 7 max = 2 x 10 ^ which results 
in high energy tail of jet emission to match the Fermi/hAT 
spectrum. We note that the choice of 7 max is very depen¬ 
dent on preferred 7 -ray emission mechanism i.e. if SSC was 
the dominant component at MeV-GeV energies reproduc¬ 
ing Fermilh AT spectra would r equire much larger value of 
7max ~ 10® (iTanaka et ahllioisl j. 

We find that choosing jet power Ljet ~ 1.3 x 10"^® erg s“® 
is sufficient to match calculated spectra with the data. Ob¬ 
tained value of jet power is smaller than accretion power 
by a factor of ^ 1 0. This is in contrast with result by 
iTanaka et al.l (l2015l l that Ljet > Lace. However, we note 
that for our calculations we used much higher value of ac¬ 
cretion disk radiative power as well as our model favours 
ERC as primary source of 7 -ray emission. 

The choice of the magnetic field strength is dictated 
by the normalisation of radio data. We set the magnetic 
field B' = 0.2 G at 0.1 pc from the central BH. Such value 
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results in magnetic field energy flux Lb ~ 4.5 x ergs“^. 
Assumed model parameters lead to kinetic energy of cold 
protons Lp ~ 9 x ergs”^ therefore jet magnetisation 
parameter o = Lb/Lp = 0.05. This result indicate that in 
the emission region at sub-parsec scales relativistic jet in 
3C 120 is already dominated by kinetic flux i.e. transition 
from initially Poynting dominated jet to matter dominated 
flow must have occurred closer to the central BH. 

Figure [3 presents broadband 3C 120 spectra for 
all Lermi/LAT data calculated with described numerical 
model. 

Since modelling “orphan” flares in 3C 120 requires an 
additional spine component we assume that the layer forms 
the underlying steady part of the jet in such a way that 
the total jet radiation consists of layer emission with ad¬ 
ditional spine emission. Note however, that during quies¬ 
cence period the spine may also exist in the jet yet due 
to large F® and large viewing angles its radiation is Doppler 
de-boosted i.e. even though the spine carries significant part 
of the total jet power it remains “invisible”. However, sud¬ 
den change of direction of the spine towards the observer 
results in sudden increase of the Doppler factor so that 
the spine contribution to total jet radiation becomes signif¬ 
icant. Simplicity of such model implies only few additional 
parameters to be determined to match observed spectra, 
namely: spine Lorentz factor F®, spine opening angle 6^ and 
an angle towards the observer ^obs- We do not make any 
further assumptions concerning the magnetic held in spine 
except for keeping spine magnetisation identical to magneti¬ 
sation of the layer. This attitude is dictated by model sim¬ 
plicity i.e. except for the spine’s magnetic held substantially 
higher than layer’s, total radiative output is dominated by 
ERC emission due to much higher Lorentz factor thus syn¬ 
chrotron and SSC components are greatly weakened in spine 
emission. 

We choose aforementioned spine parameters so that the 
calculated spectra matches the Fermi/LAT data during out¬ 
bursts 1 and 2. For hare 1 we choose = 20 and for hare 
2 we choose F®’^ = 40. We choose spine opening angles so 
that they follow a relation 9^ = O.S/F® i.e. slightly narrower 
with respect to the Lorentz factor than assumed layer open¬ 
ing angle. Such choice results in 9^’^ = 0.04 and 0”’^ = 0.02 
during hares 1 and 2, respectively. We also assume that the 
spine is bent towards the observer so that it is seen exactly 
on the border of the Doppler cone i.e. = 9^’^’^. Re¬ 

sulting Doppler factors are 22“’^ « 25 and T>^’^ « 50. Such 
choice of F® and parameters locates the emission region 
in the spine at « 1.4 x 10^^ cm « 0.04 pc - this value is 
almost identical for both hares despite signihcant diherence 
in 7 -ray hux and assumed spine Lorentz factor. Estimated 
location of the blazar zone in the spine is just outside BLR 
thus it enables avoiding signihcant absorption due to 77 pair 
production. With such parameters we are able to reproduce 
source spectra during haring states. Results of spectral en¬ 
ergy distribution (SED) modelling are presented in Fig. [S] 

Calculated spine powers are much lower than power car¬ 
ried by the layer. In case of hare 1 the ratio between spine 
and layer power Lg/Li « 0.14 and « 0.10 for hare 2 i.e. the 
layer carries majority of the jet power. 


4 CONCLUSIONS 

Fermi/LAT data analysis suggests that although it is a 
rather faint 7 -ray source it is worth monitoring especially 
as it occasionally shows very luminous hares. Detailed data 
analysis of two recently reported outbursts (see Fig. El 
proved not only extraordinary 7 -ray hux increase (by a fac¬ 
tor of ~ 50 for hare 1 and ~ 150 in case of hare 2) but also 
very short timescales of this phenomenon. As we argue in 
section 3.2 it is rather impossible to explain these factors 
with a one-zone emission model. We demonstrated that an 
introduction of the fast, narrow and wiggling spine compo¬ 
nent is a model which meets all the required properties. 

Such model proves to be very robust in explaining 
3C 120 haring behaviour as it can account for powerful 7 - 
ray hux increases without signihcant hux changes in other 
wavelengths while keeping the emission region small. Model 
robustness is also strengthened by the fact that the intro¬ 
duction of second jet component does not lead to doubling 
of model parameters. Indeed, the spine-layer scenario in¬ 
troduces only two additional model parameters: the spine 
Lorentz factor F” and its opening angle 9^. However, these 
parameters are merely just spine hux normalisation con¬ 
stants and there is a need for a viable theoretical model that 
could put constraints on their values and relations towards 
other jet properties. 

Fositive detection for time bins as short as one hour (or 
even less in case of hare 2 ) must lead to extremely small 
emitting zone if interpreted as coherent emission from lo¬ 
calised active region in a horizontally uniform jet. This ob¬ 
servation inevitably localises the 7 -ray active region very 
close to the central BH unless the jet is very fast and point¬ 
ing towards the observer which leads to high Doppler factors 
and, therefore, emission region in the jet co moving frame is 
much larger. 

We model the 3C 120 SED assuming its total radiation 
consists of the accretion-related and jet-related components. 
In case of all data, averaged Fermi/ LAT spectrum can be 
explained as having pure SSC origin llTanaka et akll^lSli or 
as being dominated by ERC radiation making use of strong 
BLR and HDR radiation in the jet environment. Note that 
neither of these models is favoured as the choice between the 
dominant 7 -ray production mechanism is very dependent 
on the Doppler factor a nd observation s indic ate that this 
value is highly variable. ICasadio et al.l (l2015l l claims that 
the observed 7 -ray flux depends not only on the energetics 
of emission region but also on the changing jet orientation 
with respect to the observer. This observation supports the 
idea of the “wiggling” jet. We note that, although such con¬ 
clusion has been made by analysing VLBI radio data, the 
same conclusion may be true in the case of the fast, narrow 
spine. In general, both the layer and the spine may change 
their directions, but due to large difference in their Lorentz 
factors observational effects are much more prominent in the 
case of the spine. 

We model 3C 120 flares 1 and 2 by varying spine’s 
Lorentz factor and by changing its direction towards the 
observer. The increase of F factor in spine by four times, 
in case of flare 1 , and by eight times, in the case of flare 
2, with respect to the average value of 5.3 are enough to 
properly reproduce the observed 7 -ray spectra. We note 
that this applies not only to flux normalisation but also 
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V [Hz] 

Figure 7. 3C 120 broadband spectrum; red points are Fermi/hAT data points for all Fermi/hAT data (from August 4th, 2008 to 
May 26th, 2015; see section 2) where arrows indicate 95% upper limits and light-grey “butterfly” plot presents best-fit power law 7 -ray 
spectra. Orange points indicate optical and X-ray data from August 2014 (see section 3.2) a nd grey poin t s are archived data collected 
from NED database. Solid blue line presents radio-loud quasar radiation template taken from lShauEf et ^ 1I2OIII I. BLR and HDR stand 
for ERG radiation on BLR and HDR seed photons, respectively. Red solid line corresponds to total jet radiative output and black solid 
line is a sum of accretion-related and jet radiation. 




V [Hz] 


u [Hz] 


Figure 8. 3C 120 broadband spectrum for flare 1 {left panel) and flare 2 {right panel); red points are Fermi/hAT data points for flares 
1 and 2 where arrows indicate 95% upper limits and light-grey “butterfly” plot presents best-fit power law 7 -ray spectra for all Fermi/LAT 
data. Orange points indicate optical and X-ray data from August 2014 (see section 3.2) a nd grey point s are a rchived data collected from 
NED database. Solid blue line presents radio-loud quasar radiation template taken from ishang et ^ Hoill L BLR and HDR stand for 
ERG radiation on BLR and HDR seed photons, respectively. Gray solid line presents layer radiation, red solid line corresponds to total 
spine radiative output and black solid line is a sum of accretion-related and jet (spine-flayer) radiation. 


reproduces the 7 -ray spectral shape. During flares 3C 120 
7 -ray spectrum was much harder (F ~ 2 . 0 ) that the av¬ 
eraged value (F ~ 2.7). The characteristic energy of BLR 
photons (~ 10 eV) and HDR photons (~ 0.6 eV) to¬ 
gether with calculated electron distribution spectral brake 
7 b and shift towards higher energies by larger F° localises the 
ERC(BLR-I-HDR) spectral peak in the right position match¬ 
ing the observations. We also point out that both in all data 


spectra and during flares 3C 120 is visible by Fermi/hAT 
only up to energies of about several GeVs. This may indi¬ 
cate that a part of the high energy radiation is absorbed via 
77 pair production on dense BLR photon field. 

The calculated powers carried by spine and layer in 
the form of kinetic energy of cold protons and Poynting 
flux indicate that majority of the total jet energy is con¬ 
tained in the layer. This is in contrast to models proposed 
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Table 2. Parameters used in numerical simulations. Single, centred values are common for all FermilhKT and flare 1 & 2 data. 


Parameter 

All data 
(layer) 

Flare 1 Flare 2 

(spine) 

Black hole mass Mbh 


4.6 X 10®A7q 


Accretion rate M 


0.3 hEdd/c^ 


Accretion disk radiative efficiency 77disk 


0.1 


Energy dissipation efficiency ?7diss 


0.3 


Jet Lorentz factor P 

5.3 

20.0 

40.0 

Doppler boosting factor X> 

3.7 

25 

50 

Fraction of energy transferred to electrons r]e 


0.5 


Magnetic field B at 0.1 pc [G] 

0.2 

0.1 

0.25 

Jet magnetisation a 


0.05 


Pair content Ue/up 


0.5 


Electron injection function indices pi,p 2 


0.5, 2.4 


Min. and max. injection energies 7mirn7max 


1, 2 X 10"^ 


Jet opening angle 

0.19 

0.04 

0.02 


l/F 

i/(o.8r) 

i/(o.8r) 

Observing angle 0obs 

0.26 

0.04 

0.02 

Location of emission region r 

2.6 X lOi'^ cm 9.5 x lO^® cm 

7.3 X lOi® cm 

BLR photons energy i^blr 


10 eV 


HDR photons energy i^'hdr 


0.06 - 0.6 eV 


BLR radius tblr 


6.5 X lO^® cm = 0.02 

pc 

HDR radius thdr 


6.5 X lO^’^ cm = 0.2 

pc 

BLR covering factor ^blr 


0.1 


HDR covering factor ^hdr 


0.3 


Accretion disk luminosity Ldisk [^rg s“^] 

Jet power Ljet [erg s“^] 

Kinetic energy of protons Lp [erg s“^] 
Magnetic energy flux Lb [erg s“^] 

1.3 X lO'^® 
9 X 10'*^ 
4.5 X lO'^® 

1.7 X 10^® 

1.8 X lO"*"* 

1.2 X lO"*"* 

6.0 X 10"^^ 

1.3 X lO"*"* 

9.0 X 10-*® 

4.5 X 10-*® 


bv iGhisellini et ^ (l2005l l where most of the total jet power 
was carried by the spine exceeding layer’s power by at least 
several times. 

Spine Lorentz factors proposed by our modelling are 
mnch different from average value of 5.3 so they are differ¬ 
ent from each other. We point out that the proposed model 
which assumes different F® factors for flares 1 and 2 is not 
unique i.e. it is possible to get similar results with differ¬ 
ent set of input parameters. It is especially important to 
notice that the observed radiation is very sensitive to as¬ 
sumed Doppler factors, hence, the direction of the fast spine 
towards the observer. It is very plausible that both flares 
could be modelled with identical spine Lorentz factor with 
different direction towards the observer as the main reason 
for different observed 7 -ray fluxes. 

Another way of having such extreme values in a sin¬ 
gle source is that the spine-layer model might be just a 
coarse simplification and yet an indication of a much more 
complex jet structure. Instead of a sharp division into fast 
and slow components one may think of a smooth transition 
from slower, outer parts of the jet towards faster, inner re¬ 
gions e.g. described by a Gaussian profile r(i?) where R is 
the distance from the jet axis. Depending on the broadness 
of such Gaussian F factors distribution observed radiation 
could be characterised by different Lerc/Lsyu ratio. Also 
the variability of this ratio would be very sensitive on the jet 
F(i?) profile i.e. small short-term 7 -ray variability would in¬ 
dicate rather broad prohle while rapid, powerful high energy 
flares could be explained with peaked, narrow distribution 
similar to model proposed in this work. However, detailed 


and quantitative model of such proposition is required to 
assess its validity. 
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